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This  report  presents  the  results  of  a  study  conducted  for 
the  Federal -State  Interagency  San  Joaquin  Valley  Drainage 
Program.  The  purpose  of  the  report  is  to  provide  the  Drainage 
Program  agencies  with  information  for  consideration  in 
developing  alternatives  for  agricultural  drainage  water 
management.  Publication  of  any  findings  or  recommendations  in 
this  report  should  not  be  construed  as  representing  the 
concurrence  of  the  Program  agencies.  Also,  mention  of  trade 
names  or  commercial  products  does  not  constitute  agency 
endorsement  or  recommendation. 


The  San  Joaquin  Valley  Drainage  Program  was  established  in 
mid-1984  as  a  cooperative  effort  of  the  U.S.  Bureau  of  Reclamation, 
U.S.  Fish  and  Wildlife  Service,  U.S.  Geological  Survey,  California 
Department  of  Fish  and  Game,  and  California  Department  of  Water 
Resources.  The  purposes  of  the  Program  are  to  investigate  the 
problems  associated  with  the  drainage  of  irrigated  agricultural  lands 
in  the  San  Joaquin  Valley  and  to  formulate,  evaluate,  and  recommend 
alternatives  for  the  immediate  and  long-term  management  of  those 
problems.  Consistent  with  these  purposes,  Program  objectives  address 
the  following  key  areas:  (1)  Public  health,  (2)  surface-  and  ground- 
water resources,  (3)  agricultural  productivity,  and  (4)  fish  and 
wildlife  resources. 

Inquiries  concerning  the  San  Joaquin  Valley  Drainage  Program  may 
be  directed  to: 


San  Joaquin  Valley  Drainage  Program 
2800  Cottage  Way,   Room  W-2143 
Sacramento,  California    95825-1898 
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EXECUTIVE  SUMMARY 

PHASE  II  REPORT  TO  THE  SAN  JOAQUIN  VALLEY  DRAINAGE  PROGRAM: 

STUDY  OF  INNOVATIVE  TECHNIQUES  FOR  SUBSURFACE  DRAINAGE  REDUCTION 


This  report  summarizes  results  of  the  first  year  of  field  studies 
of  irrigation  management/controlled  drainage  techniques  to  reduce 
subsurface  drainage  flows  and  constituent  loads.  Three  furrow- 
irrigated  fields  with  tile  drains  were  selected  in  the  Grasslands 
Basin  for  the  study.  Two  fields  were  in  the  Panoche  Water  and 
Drainage  Districts.   The  third  was  in  Broadview  Water  District. 

The  first  field  in  Panoche  Water  District,  Field  P-1,  was  a  93  acre 
cotton  field.  During  the  1988  crop  season  an  irrigation  management 
program  was  implemented.  Tile  drain  weir  valves  were  installed, 
and  set  to  control  the  drainage  outflow  from  the  field.  Irrigation 
and  drainage  monitoring  data  showed  that  the  combination  of 
irrigation  management  and  controlled  drainage  resulted  in  limiting 
net  deep  percolation  to  3.5  inches.  This  was  approximately  equal 
to  the  minimum  net  deep  percolation  amount  which  must  be  achieved 
to  maintain  the  field's  salt  balance.  The  tile 'drain  outflow  also 
amounted  to  3.5  inches.  Thus,  over  the  course  of  the  season, 
lateral  and  deep  subsurface  inflows  balanced  lateral  and  deep 
subsurface  outflows.  The  tile  drain  outflow  from  individual  tile 
drains  varied  significantly.  This  could  be  attributed  to 
variations  in  the  tile  drain  depths.  Salt  loads  into  the  field 
from  irrigation  water  were  roughly  equal  to  salt  loads  leaving  the 
field  in  the  tile  drain  water.  Selenium  levels  in  the  tile  drain 
water  were  50%  higher  than  those  in  the  irrigation  water.  Boron 
loads  from  the  tile  drains  were  more  than  three  times  the  amount 
which  infiltrated  with  the  irrigation  water. 

The  second  field  in  Panoche  Water  District,  Field  P-2,  was  a  150 
acre  cotton  field.  During  the  1988  crop  season  an  irrigation 
management  and  drainage  monitoring  program  was  initiated.  A  lack 
of  monitoring  history,  large  water  table  variations  and  steep  field 
slope  prevented  controlled  drainage  from  being  implemented  in  1988. 
The  goal  of  this  year's  program  was  to  gain  sufficient  information 
about  the  field's  dynamics  to  allow  controlled  drainage  in  1989. 
Irrigation  and  drainage  monitoring  data  showed  that  even  with  the 
irrigation  management  efforts,  the  net  deep  percolation  amounted 
to  18  inches.  This  exceeded  the  minimum  deep  percolation  by  16 
inches.  This  was  attributed  mostly  to  the  high  intake  nature  of 
the  field's  soils.  Eleven  inches  of  water  flowed  out  of  the 
subsurface  drains.  The  tile  drain  outflows  from  the  individual 
tile  drains  varied  significantly.  This  could  be  attributed  mostly 
to  the  variations  in  water  table  levels  observed  in  the  field. 
Seven  inches  of  water  were  lost  in  deep  and  lateral  subsurface 
flows.  Tile  drain  salt,  boron  and  selenium  loads  were  1.5,  5.1  and 
1.3  times  the  amount  of  salt,  boron  and  selenium  infiltrating  with 
the  irrigation  water. 


The  field  in  Broadview  Water  District,  Field  B10-2,  was  a  150  acre 
melon  field.  During  the  1988  crop  season,  an  irrigation  management 
and  drainage  monitoring  program  was  initiated.  Control  of  the  tile 
drainage  system  with  weir  valves  was  postponed  until  1989  in  order 
to  develop  a  better  understanding  of  the  field's  dynamics. 
Irrigation  and  drainage  monitoring  data  showed  that  the  irrigation 
management  efforts  resulted  in  limiting  net  deep  percolation  to 
four  inches.  This  was  still  more  than  the  minimum  net  deep 
percolation.  Tile  drain  outflows  amounted  to  4.8  inches.  Incoming 
deep  and  lateral  subsurface  flows,  therefore,  exceed  outgoing  deep 
and  lateral  subsurface  flows  by  approximately  one  inch.  Tile  drain 
salt,  boron  and  selenium  loads  were  3.7,  3.8  and  3.9  times  the 
amount  of  salt,  boron  and  selenium  infiltrating  with  the  irrigation 
water  respectively. 

DRAINMOD  simulations  based  upon  field-gathered  input  data  were  run 
for  comparison  to  actual  water  table  and  tile  drain  flow  data  for 
each  field.  A  general  correlation  between  DRAINMOD  simulations  and 
actual  field  data  was  observed  for  all  three  fields.  Significant 
differences  between  DRAINMOD  and  actual  data  could  be  attributed 
mostly  to  influences  from  neighboring  fields. 

The  general  correlation  between  DRAINMDD  and  actual  data  suggests 
that  DRAINMOD  will  be  useful  for  analyzing  the  management  impacts 
of  various  field  activities.  DRAINMOD's  ability  to  predict  actual 
field  conditions  currently  depends  upon  the  strength  of  outside 
influences.  A  newly  released  version  of  DRAINMOD  may  help  model 
some  of  the  outside  influences. 

Recommendations  are  given  on  modifications  which  are  needed  to 
DRAINMOD  to  make  it  easier  to  use  and  more  applicable  to  the 
management  of  irrigated  agriculture. 


A  further  refinement  to  the  controlled  drainage  program  to  be 
implemented  next  year  will  be  to  set  the  weirs  based  upon  depth 
below  the  ground  surface  rather  than  height  above  tile  inverts. 

Results  of  the  field  studies  to  date  have  strengthened  JMLord's 
opinion  that  irrigation  management/controlled  drainage  techniques 
can  be  used  to  reduce  subsurface  drainage  flows  and  loads. 
Controlled  drainage  in  the  San  Joaquin  Valley  is  an  emerging 
technology  and  more  work  is  needed  to  develop  it  into  a  service 
which  can  be  applied  on  a  large  scale.  Preliminary  estimates  of 
hardware  and  management  costs  for  irrigation  management/controlled 
drainage  are  given. 
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SJVDP  PHASE  II  REPORT 
STUDY  OF  INNOVATIVE  TECHNIQUES  FOR  DRAINAGE  REDUCTION 


I.   INTRODUCTION 

A.  Study  goals 

This  report  summarizes  results  from  the  first  year  of  field 
investigations  of  irrigation  management  and  controlled  drainage 
techniques.  The  purpose  of  these  investigations  is  to  reduce 
subsurface  drainage  flows  and  their  constituent  loads. 

This  report  covers  the  second  phase  of  a  three  year  study  for  the 
San  Joaquin  Valley  Drainage  Program  (SJVDP).  The  first  phase  was 
to  study,  in  theory,  the  concept  of  irrigation  management  with 
controlled  drainage.  This  phase  was  completed  with  a  report 
printed  in  November  1987  [1].  The  second  and  third  phases  of  the 
study  were  to  develop  the  management  techniques  and  demonstrate 
their  implementation  on  three  fields  with  existing  tile  drainage. 
All  of  these  fields  are  in  the  Grasslands  Basin. 

The  fields  were  located  and  the  cooperation  of  the  growers  was 
obtained  in  late  1987.  Figure  1  is  a  map  showing  the  location  of 
the  three  fields.  Table  1  presents  a  summary  of  the  field 
information.  A  water  table  control  strategy  was  implemented  on 
one  of  the  fields  with  the  other  two  fields  being  put  under  a  water 
management  and  monitoring  program.  A  great  deal  of  progress  was 
made  toward  being  able  to  implement  water  table  control  in  all 
three  of  the  fields  in  the  coming  year  for  phase  three. 

B.  Field  P-1 

Field  P-1  has  a  net  area  of  91.7  acres.  It  was  planted  to  cotton 
in  1988  following  a  wheat  crop  in  1987.  It  is  located  within  the 
Panoche  Water  and  Drainage  Districts.  Figure  2  is  a  layout  of  the 
field  with  the  associated  irrigation  and  drainage  facilities. 

This  field  has  six  buried  plastic  tile  drain  laterals  (installed 
in  1978)  which  are  part  of  a  much  larger  network  of  tile  drains 
for  several  fields  in  the  surrounding  area.  At  the  beginning  of 
the  study,  the  farm  provided  a  map  which  showed  all  of  the  tile 
drains  being  parallel  and  evenly  spaced  400  feet  apart.  After  the 
field  was  selected  for  the  study,  an  old  aerial  photo  was  examined. 
It  was  discovered  (to  the  farm's  surprise)  that  the  drains  were  not 
installed  parallel  and  were  not  evenly  spaced.  Furthermore,  when 
tile  drain  valves  were  installed  in  each  tile  line,  it  was 
discovered  that  three  of  the  drains  were  installed  about  7  feet 
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deep  and  three  were  about  8  feet  deep.  All  of  the  drains  had  been 
stretched  during  installation,  were  oval  shaped,  and  had 
significant  amounts  of  sediment  in  the  pipes.  Consideration  was 
given  to  looking  for  another  field  with  more  uniform  conditions  to 
study.  However,  the  decision  was  made  to  stay  with  this  field 
because  the  non-uniform  tile  system  is  likely  typical  of  many  in 
the  area.  The  conclusion  made  was  that  if  controlled  drainage 
worked  in  this  field,  it  would  be  reasonable  to  assume  that  the 
concept  would  work  in  other  fields,  both  non-uniform  and  uniform. 

There  is  a  4.5  foot  fall  in  the  ground  surface  elevation  from  South 
to  North.  There  was  a  question  of  whether  one  valve  per  tile  line 
would  be  enough  to  control  a  field  with  that  much  fall.  It  was 
decided  to  try  one  valve  in  1988  and  go  to  two  valves  in  1989  if 
one  was  determined  to  not  be  enough.  The  tile  drain  valves  were 
installed  at  the  north  edge  of  the  field.  These  valves  are 
inserted  in  the  tile  line  and  are  configured  such  that  two  movable 
plates  function  as  weirs  to  control  the  upstream  water  level.  The 
weirs  were  set  to  30  inches  above  the  drain  invert  just  prior  to 
pre-irr igation.  After  the  second  in-season  irrigation,  the  weirs 
were  set  to  the  ground  surface  to  totally  resist  all  flow. 

C.   Field  P-2 

Field  P-2  has  a  net  area  of  149.9  acres.  It  was  planted  to  cotton 
in  1988  following  a  safflower  crop  in  1987.  The  field  is  within 
the  Panoche  Water  and  Drainage  Districts.  Figure  3  shows  the 
irrigation  and  drainage  features. 

There  are  six  buried  concrete  drain  tile  laterals  spaced  400  feet 
apart  and  6  feet  deep  which  run  the  length  of  the  field  and 
discharge  into  a  deep  open  drain.  The  field  is  irrigated  in  one 
quarter  mile  runs.  The  west  half  is  irrigated  by  a  head  ditch  with 
siphons.  The  east  half  is  irrigated  by  a  line  of  gated  pipe.  The 
ground  surface  elevation  slopes  from  the  Southwest  to  the 
Northeast.  From  the  west  edge  to  the  east  edge  there  is  12  feet 
of  fall.  Because  of  this,  controlling  the  drain  flows  would 
require  a  valve  every  660  feet  (l/8th  mile)  in  order  to  have  3  feet 
of  fall  from  valve  to  valve.  However,  the  irrigation  deliveries, 
drainage  outflows,  and  water  quality  of  the  field  had  never  been 
monitored,  so  there  was  no  way  to  compare  the  results  of  a 
controlled  drainage  trial  with  previous  years.  It  was  decided  to 
install  one  tile  drain  valve  in  each  of  the  tile  lines  at  mid-field 
prior  to  the  1988  pre-irr igation.  The  weirs  would  not  be  set  to 
control  the  outflows  but  would  be  used  as  an  access  point  for  flow 
measurement.  The  1988  season  would  be  used  to  start  work  on  the 
irrigation  management  aspect  of  the  program,  for  gathering  data  on 
the  subsurface  drainage  system  for  later  comparison,  and  for 
developing  a  management  strategy  for  the  1989  crop  season  using 
controlled  drainage. 
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D.   Field  B10-2 

Field  B10-2  has  a  net  area  of  145.5  acres.  It  was  planted  to 
cantaloupes  in  1988  following  a  cotton  crop  in  1987.  The  field  is 
within  the  Broadview  Water  District.  Figure  4  shows  the  irrigation 
and  drainage  features. 

This  field  has  nine  buried  corrugated  plastic  tile  laterals 
connected  to  a  concrete  tile  baseline.  The  field  to  the  West 
(Field  B10-1)  also  has  nine  tile  laterals  which  flow  into  the  same 
base  collector  line.  The  tile  base  line  ends  near  the  northeast 
corner  of  the  field  at  a  steel  sump  with  a  pump.  The  pump 
discharges  the  tile  drain  water  into  a  Broadview  Water  District 
open  drain  that  runs  along  the  east  edge  of  the  field.  The  outflow 
from  the  sump  has  been  monitored  by  the  Broadview  Water  District 
for  flow,  selenium  concentration,  and  dissolved  salts  since  1985. 

This  field  is  irrigated  in  one  quarter  mile  runs.  Both  halves  are 
irrigated  by  a  head  ditch  with  siphons.  The  ground  surface 
elevation  slopes  from  the  Southwest  to  the  Northeast.  From  the 
south  edge  to  the  north  edge  there  is  5  feet  of  fall.  To 
adequately  control  the  water  level  with  this  slope  necessitates 
installing  a  valve  in  each  tile  line  at  the  midpoint. 

Since  melons  are  relatively  sensitive  to  salt,  boron,  and  saturated 
conditions,  JML  felt  that  it  was  too  risky  to  try  to  implement  a 
severely  managed  drainage  program  without  first  getting  a  better 
understanding  of  the  field's  subsurface  drainage.  For  this  reason, 
it  was  decided  that  the  1988  season  would  be  used  to  start  work  on 
the  irrigation  management  aspect  of  the  program  and  to  develop  a 
better  understanding  of  the  dynamics  of  the  field.  The  goal  would 
be  to  develop  a  management  strategy  using  controlled  drainage  for 
the  1989  crop  season.  In  order  to  measure  the  flow  of  water  from 
Field  B10-1  separately  from  B10-2,  a  tile  drain  valve  was  installed 
in  the  mainline  as  an  access  hole  for  a  flow  meter. 
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II.   METHODS 


A.  Irrigation  deliveries  and  tailwater  flows 

Irrigation  deliveries  to  each  of  the  fields  were  measured  using 
propeller  flow  meters  with  instantaneous  and  totalizing  indicators. 
Water  quality  tests  of  the  irrigation  water  were  done  by  the 
districts. 

The  tailwater  runoff  from  Fields  P-1  and  B10-2  was  measured  with 
Parshall  flumes  with  continuous  water  level  recorders  in  the 
stilling  wells.  Following  each  irrigation,  a  tailwater  hydrograph 
was  developed.  At  Field  P-2,  nearly  all  of  the  tailwater  is 
recycled  for  reuse  on  the  same  field.  There  is  a  spill  pipe  in  the 
tailwater  pit  which  spills  water  into  the  T  Canal  Drain  during  peak 
tailwater  flows.  This  spill  was  monitored  by  installing  a  stilling 
well  with  a  continuous  water  level  recorder  on  the  tailwater  pit 
side  of  the  spill  pipe  .  Spot  checks  of  the  spill  flow  rate  were 
made  with  a  bucket  and  stop  watch.  From  this  a  water  level  vs. 
spill  flow  rate  relationship  was  developed.  This  relationship  was 
used  to  convert  the  recorder  charts  into  hydrographs. 

Evaluations  were  done  during  each  of  the'  first  three  irrigations 
on  each  field  to  estimate  the  uniformity  of  distribution  of 
infiltrated  irrigation  water.  A  two  point  volume  balance  method 
was  used.  This  method  was  based  upon  the  furrow  irrigation 
theories  of  Walker  and  Skogerboe  of  Utah  State  University  [2]. 
The  same  theory  is  the  basis  of  the  "Waterman  Program"  used  in  the 
Westside  Resource  Conservation  District  and  Westlands  Water 
District's  Water  Conservation  and  Drainage  Reduction  Program. 

B.  Tile  drain  flows  and  quality 

Three  different  measurement  techniques  were  used  to  measure  the 
flow  of  tile  drain  water.  At  Field  P-2,  the  gravity  outfalls  to 
the  deep  drain  were  measured  using  a  graduated  five  gallon  bucket 
and  a  stop  watch.  This  method  was  very  accurate  for  flows  in  the 
1  to  60  gpm  flow  range.  Higher  flows,  which  occurred  for  short 
periods  of  time  following  irrigations,  were  less  accurate. 
Therefore,  readings  were  repeated  several  times  and  averaged  during 
high  flows  to  reduce  the  error. 

At  all  three  fields,  flows  were  measured  at  the  tile  drain  valves 
with  an  electro-magnetic  velocity  meter  (Marsh-McBirney  Model  201- 
D  current  meter).  Before  inserting  the  meter,  the  weir  was  set  to 
cause  the  inlet  pipe  to  flow  full.  After  allowing  the  water  level 
in  the  valve  to  stabilize,  the  velocity  meter  was  inserted  at  the 
mid  point  of  the  inlet  pipe  and  then  the  velocity  reading  taken. 
Velocity  readings  were  converted  to  flow  rates  by  multiplying 
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velocity  times  cross  sectional  area  times  a  conversion  factor.  The 
conversion  factor  was  developed  by  putting  an  upturned  PVC  elbow 
on  the  discharge  of  one  of  the  tile  drain  outfalls  in  Field  P-2  and 
by  taking  simultaneous  bucket  catches  and  velocity  meter  readings 
at  different  times  of  the  year. 

At  Field  B10-2,  a  propeller  meter,  which  was  installed  on  the 
discharge  of  the  sump  pump  by  Broadview  Water  District,  was  used 
to  obtain  flows  from  both  Fields  B10-1  and  B10-2. 

Since  the  bucket  and  velocity  meter  readings  gave  instantaneous 
flows  but  not  volumes,  flow  volumes  were  calculated  by  averaging 
the  flow  from  two  readings  and  multiplying  by  the  interval  between 
readings.  The  accuracy  of  the  flow  volumes  depends  upon  reading 
flows  frequently  enough  so  as  not  to  miss  any  peaks  or  valleys  in 
the  hydrograph.  The  frequency  of  the  flow  readings  taken  during 
the  study  was  designed  to  be  flexible  according  to  the  amount  of 
change  which  was  taking  place.  During  and  shortly  after 
irrigations,  readings  were  taken  every  two  to  three  days. 
Following  this,  the  reading  interval  was  gradually  increased  to  a 
maximum  of  two  weeks.  This  flexible  schedule  allowed  good  data  to 
be  gathered  while  minimizing  monitoring  man-hours. 

Monthly  grab  samples  were  taken  from  each  measuring  point  for 
laboratory  analysis.  Each  sample  was  tested  for  electrical 
conductivity  (EC)  and  boron  in  JMLord,  Inc.'s  laboratory.  EC 
analysis  was  done  with  a  Wheatstone  bridge  meter  with  automatic 
temperature  adjustment.  Boron  analysis  was  done  using  the 
colorimetric  method.  Each  sample  was  also  tested  for  selenium 
using  the  hydride  generation/atomic  absorption  spectrophotometry 
method.  Some  of  the  selenium  samples  were  run  by  Twining 
Laboratories  in  Fresno.  Others  were  run  by  the  USDA's  Water 
Management  Laboratory  in  Fresno. 

C.   Shallow  groundwater  levels  and  quality 

The  shallow  groundwater  levels  in  each  field  were  monitored  by 
installing  10  foot  long,  slotted,  1-1/2  inch  diameter  PVC  pipe  in 
9-1/2  foot  deep  auger  holes.  Each  hole  was  backfilled  with  San 
Joaquin  River  sand.  The  holes  were  drilled  with  a  3  inch  hand 
auger.  The  slots  were  cut  with  a  carpenters  saw  every  4  inches 
from  three  to  ten  feet  of  depth.  The  bottom  of  the  pipe  was  capped 
to  prevent  soil  intrusion.  The  locations  of  the  observation  wells 
were  selected  to  provide  as  much  information  about  the  spatial 
variability  of  the  groundwater  in  the  field  as  possible,  while 
limiting  the  number  of  sites  to  around  25. 

After  the  observation  wells  were  installed,  each  field  was  surveyed 
with  a  theodolite/EDM  survey  instrument.  This  was  done  to 
precisely  map  the  observation  well  locations  and  elevations.   The 

JMLord,  Inc.  Page 

Scientists  and  Engineers  6 


depth  to  groundwater  levels  were  read  on  the  same  schedule  as  the 
tile  drain  flows  with  a  simple  device  JML  calls  a  "plopper".  It 
consists  of  a  hollow  tipped  weight  on  a  marked  string.  To  read  the 
depth  to  groundwater,  a  technician  lowers  the  weight  into  the  pipe 
until  a  splash  is  heard.  The  top  of  the  pipe  is  used  as  a 
reference  and  the  markings  on  the  string  indicate  the  depth 
reading.  The  distance  from  the  top  of  the  pipe  to  the  ground 
surface  level  (obtained  in  the  survey)  is  then  subtracted  from  the 
reading  to  obtain  the  true  depth  to  groundwater. 

The  observation  wells  were  also  used  as  points  for  withdrawing 
shallow  groundwater  for  laboratory  analysis.  Samples  were  taken 
once  in  the  winter  and  once  in  the  summer  from  all  of  the 
observation  wells  and  analyzed  for  EC  and  boron.  Samples  were 
obtained  by  dropping  a  weighted,  elongated  cup  into  the  wells  and 
withdrawing  the  sample. 

D.  Soil  physical  and  chemical  properties 

Each  hole  was  logged,  as  it  was  being  augured,  for  soil  texture, 
structure,  apparent  permeability,  color,  and  moisture  content.  In 
addition,  samples  from  each  foot  from  the  top  five  feet  of  each 
well  site  were  put  into  separate  bags  for  laboratory  analysis.  The 
pH ,  EC,  boron,  sodium,  calcium,  magnesium,  bicarbonate,  carbonate, 
chloride,  and  nitrate  levels  were  analyzed  on  saturated  soil  pastes 
in  JML ' s  laboratory.  In  addition,  one  time  during  the  summer,  a 
diagonal  pass  was  made  across  each  field  to  collect  composite 
samples  from  each  of  the  top  three  feet  of  soil  for  the  same  set 
of  analyses. 

Based  upon  the  soil  profile  logs  made  while  auguring,  four  sites 
were  chosen  in  each  field  for  auger  hole  hydraulic  conductivity 
measurements.  These  were  based  on  the  bailout  method  described  in 
the  Bureau  of  Reclamation's  Drainage  Manual  [3].  The  tests  were 
run  in  four-inch  auger  holes  shortly  following  the  pre-irr igation 
of  each  field. 

Core  samples  from  each  foot  of  the  top  five  feet  of  soil  from  the 
hydraulic  conductivity  test  holes  were  taken  for  laboratory 
analysis  of  the  soil  moisture  characteristic  curves.  The  soil 
moisture  characteristic  curves  were  obtained  in  JML '  s  laboratory 
by  measuring  the  moisture  content  of  each  soil  while  saturated  and 
at  pressures  of  0.1,  0.333,  1,  5,  and  10  bars  on  a  pressure  plate 
apparatus. 

E.  Evapotranspirat ion  and  soil  moisture. 

Evapotranspiration  rates  were  calculated  using  weather  data  from 
the  Telles  weather  station  operated  by  the  California  Irrigation 
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Management  Information  Service  (CIMIS).  A  modified  Penman  equation 
was  used  by  JML ' s  in-house  computer  irrigation  scheduling  program 
(SWAP/ET)  to  calculate  the  evapotranspirat ion  potential  (ETP) . 
SWAP/ET  multiplies  ETP  by  crop  factors  taken  from  crop  curves  to 
obtain  estimated  evapotranspirat ion. 

During  the  crop  season,  soil  moisture  levels  were  measured  with  a 
neutron  probe  at  three  different  sites.  Readings  were  taken  from 
zero  to  five  feet,  in  one-foot  increments,  at  each  site  before  and 
after  each  irrigation  and  every  week  between  irrigations.  The 
probe  was  calibrated  using  gravimetric  samples  taken  from  the 
access  tube  sites  as  they  were  installed. 


F.   DRAINMOD  simulations 

In  order  to  manage  the  irrigations  and  control  the  weir  settings 
as  JML  plans  to  do,  it  is  important  to  be  able  to  predict  what 
would  happen  to  the  water  table,  drain  flows,  soil  moisture,  etc. 
in  a  tile  drained  field  under  different  management  strategies. 
The  tool  that  JML  chose  for  this  purpose  is  DRAINMOD.  DRAINMOD  is 
a  computer  program  which  models  humid  zone  drainage  hydrology  [4]. 
It  has  not  been  extensively  tested  in  irrigated  arid  regions. 
During  the  1988  crop  season  the  field  data  necessary  to  allow  the 
testing  of  DRAINMOD  as  a  management  tool  was  gathered.  The  goal 
was  to  be  able  to  make  the  predictions  which  are  necessary  for 
determining  the  future  irrigation  management  and  controlled 
drainage  techniques  to  be  used  in  the  three  fields. 

The  drainage  system  configurations,  initial  water  table  depths, 
soil  moisture  characteristic  curves,  hydraulic  conductivity  test 
results,  actual  irrigation  schedules,  rainfall  events,  and  the 
estimated  daily  ET  values  from  the  data  collected  were  used  as 
inputs  to  the  model  for  testing.  Results  of  the  DRAINMOD 
simulations  were  then  compared  to  the  collected  field  data  for 
depths  to  water  table,  tile  drain  flows,  soil  moisture  content, 
and  upflux. 


JMLord,  Inc.  Page 

Scientists  and  Engineers  8 


III.   RESULTS 

A.   Field  P-1 

i.   Irrigation 

Table  2  summarizes  data  collected  on  the  irrigations  of  Field  P-1 
in  1988. 

Table  2.   Field  P-1  Irrigation  Summary 


Start 

End 

Gross 

Run- 

Avg. 

Dist . 

Date 

Date 

Applied 

off 

Infilt. 

Unif .* 

(dd/mm) 

(dd/mm) 

(in.) 

(in.) 

(in.) 

(%) 

02/09 

02/27 

15.4 

5.3 

10.1 

97 

06/09 

06/14 

4.3 

1.0 

3.3 

85 

07/04 

07/11 

5.5 

2.2 

3.3 

90 

07/22 

07/28 

4.7 

1.8 

2.9 

90# 

08/09 

08/14 

4.8 

2.4 

2.4 

90# 

Total 

for  year 

34.7 

12.7 

22.0 

92 

*  Average  distribution  uniformity  of  test  furrows 

#  Estimated  values,  as  they  were  not  measured 

The  pre-season  irrigation  was  made  in  February.  As  a  result  of 
JML ' s  evaluation  of  the  performance  on  the  first  set,  the  farmer 
changed  from  48  hour  sets  with  very  low  furrow  flows  to  24  hour 
sets  with  higher  flows  for  the  remainder  of  the  irrigation.  This 
reduced  the  amount  of  infiltrated  water  by  about  1  inch.  The 
amount  which  infiltrated  was  still  much  more  than  was  necessary  to 
replenish  the  potential  root  zone  for  the  coming  season.  However, 
the  soil  intake  characteristics  and  the  need  to  adequately  "sub" 
up  the  furrows  to  provide  seed  bed  moisture  prevented  any  further 
reductions . 

The  first  in-season  irrigation  was  done  in  June.  This  was  done  on 
every  other  furrow  to  reduce  the  amount  infiltrated.  During  the 
last  three  irrigations,  the  soil  had  sealed  up  to  where  the  soil 
moisture  was  not  being  completely  replenished. 

ii.   Tile  drain  flows,  quality,  and  loads 
Table  3  is  a  summary  of  the  tile  drain  data  for  the  1988  crop  year. 
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Table  3.   Field  P-1  Tile  Drain  Data  Summary 


Tile 

Total 

Avg . 

Avg  . 

Avg  . 

Salt 

Boron 

Se 

Number 

Flow 

EC 

Boron 

Se 

Load 

Load 

Load 

(#) 

(ac-ft) 

(dS/m) 

(mg/1) 

(mg/1) 

(tons) 

(lbs) 

(lbs) 

1 

5.9 

4.16 

7.14 

0.015 

20.9 

82.2 

0.20 

2 

0.9 

5.42 

3.67 

0.021 

4.5 

9.1 

0.05 

3 

1.2 

4.77 

4.25 

0.018 

5.2 

13.8 

0.06 

4 

3.0 

3.41 

4.50 

0.012 

9.7 

36.8 

0.10 

5 

11.5 

3.67 

7.11 

0.015 

40.3 

212.6 

0.45 

6 

5.0 

4.00 

7.41 

0.015 

19.7 

87.4 

0.20 

Total 

27.5 
(3.5  in. ) 

3.83 

5.91 

0.014 

100.2 

441.7 

1.06 

The  flow  of  tile  water  and  the  constituent  loads  are  considerably 
higher  for  the  three  deeper  tile  drains  (Numbers  1,  5  and  6) . 
Figure  5  is  a  graph  of  tile  flows  versus  time.   The  three  deeper 
tile  lines  flowed  for  a  much  longer  period  of  time  than  the  three 
shallower  ones. 

Average  laboratory  results  for  Panoche  Water  District's  water 
supply  for  1988  were:  EC  =  0.68  dS/m  (435  mg/1),  boron  =  0.31  mg/1, 
selenium  =  0.0015  mg/1.  When  these  concentrations  are  multiplied 
by  the  amount  of  infiltrated  irrigation  water,  the  resulting  loads 
are:   salt  =  100  tons,  boron  =  141  lbs,  selenium  =  0.69  lbs. 

Comparing  the  incoming  to  outgoing  loads  shows  that  the  incoming 
salts  were  equal  to  the  outgoing  salts,  outgoing  boron  was  roughly 
3  times  higher  than  incoming  boron,  and  outgoing  selenium  was 
roughly  50  percent  higher  than  incoming  selenium. 

ill.   Shallow  groundwater  levels  and  quality 

The  spatial  variation  of  the  water  table  is  illustrated  by  Figure 
6.  This  figure  shows  the  depth  to  groundwater  map  for  April  29, 
1988.  On  that  day,  the  depth  to  groundwater  varied  from  a  low  of 
5.4  to  a  high  of  7.3  feet.  Comparing  the  depth  to  groundwater  map 
with  a  ground  surface  elevation  contour  map  (Figure  7)  shows  the 
shallowest  groundwater  occurring  in  the  southeast  corner  of  the 
field,  which  is  a  low  elevation  area.  The  effect  of  the  water 
backing  up  in  the  tile  drains  behind  the  valves  did  not  appear  to 
cause  large  localized  groundwater  mounds,  although  the  water  table 
was  about  one  half  foot  shallower  on  the  north  side  of  the  field. 

The  time  variation  of  the  water  table  is  illustrated  in  Figure  8. 
This  figure  shows  a  plot  of  selected  observation  wells  through  the 
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1988  crop  season.  The  plot  shows  the  water  table  responded 
significantly  to  the  pre-season  irrigation  and  first  in-season 
irrigation  but  had  little  response  to  the  remaining  irrigations. 
The  water  table  also  showed  a  small  upward  response  during  early 
season  irrigations  on  the  tomato  fields  to  the  South. 

In  January  the  average  EC  was  4.41  ds/m  and  the  average  boron  was 
4.29  mg/1.  In  July  the  average  EC  was  6.58  dS/m  and  the  average 
boron  was  5.41  mg/1.  Figure  9  shows  the  spatial  variability  of 
the  groundwater  EC  and  Figure  10  shows  the  spatial  variability  of 
the  groundwater  boron  content.  The  mid-field  transect  of 
observation  wells  produces  a  pattern  which  shows  a  lower  quality 
water  midway  between  drains  and  a  higher  quality  water  closer  to 
the  drains. 

iv.   Soils 

Soil  profile  logs  from  the  observation  wells  indicated  that  the 
soils  in  the  top  ten  feet  of  the  field  varied  considerably.  At 
each  hole  the  soil  profile  consisted  of  layers  of  alluvial 
deposits.  However,  the  layers  did  not  appear  to  be  continuous  over 
large  areas.  The  dominant  soil  type  was  silty  clay-loam.  The  soil 
types  observed  in  the  layers  ranged  from  loamy  sand  to  dense  clay 
and  every  type  in  between.  Average  results  of  the  soil  salinity 
analysis  are  in  Table  4. 

Table  4.   Average  soil  EC  and  boron,  Field  P-1 


Depth        Avg . 
Range       ECe 
(ft)       (dS/m) 


Avg  . 

Boron 

(mg/1) 

0.94 

0.98 

1.96 

3.14 

3.68 

0-1  1.25 

1-2  1.00 

2-3  1.75 

3-4  2.55 

4-5  3.21 

Avg.  1.95        2.14 


The  soil  ECe  and  boron  varied  considerably  from  the  above  averages 
with  location.  None  of  the  sites  tested  had  average  ECe  or  boron 
levels  which  would  cause  yield  reductions  on  cotton.  Figure  11 
shows  lines  of  equal  average  root  zone  ECe  for  Field  P-1.  Figure 
12  shows  lines  of  equal  average  root  zone  boron  content  for  Field 
P-1. 

Figure  13  is  a  plot  based  upon  the  averages  of  the  soil  moisture 
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characteristic  curves  derived  for  samples  from  observation  wells 
12,  15,  19  and  26. 

In-field  tests  at  observation  wells  15,  19  and  26  gave  saturated 
horizontal  hydraulic  conductivity  values  of  2.3,  2.4,  and  4.4  in/hr 
respectively. 

V.   Evapotranspirat ion  and  soil  moisture 

Figure  14  is  a  plot  of  the  evapotranspiration  potential  (ETP)  ,  the 
ETP  multiplied  by  the  crop  coefficient  (ETP  x  Kc)  ,  and  the  crop 
evapotranspiration  (CROP  ET) .  CROP  ET  was  calculated  by  SWAP/ET 
for  scheduling  the  irrigations  and  for  use  in  DRAINMOD  simulations. 
CROP  ET  is  different  from  ETP  x  Kc.  This  is  because  a  soil  dryness 
coefficient  is  used  in  the  calculation  to  reduce  ET  when  the  root 
zone  soil  moisture  is  depleted  to  the  point  of  limiting  ET.  The 
total  growing  season  CROP  ET  was  20.0  inches.  The  total  growing 
season  ETP  x  Kc  was  25.8  inches.  Since  the  SWAP/ET  model  does  not 
include  the  effect  of  upflux  from  a  shallow  water  table,  the  actual 
soil  dryness  was  less  than  predicted  by  SWAP/ET.  For  this  reason, 
actual  ET  numbers  should  fall  somewhere  between  CROP  ET  and  ETP  x 
Kc.   For  DRAINMOD  simulations,  an  average  of  the  two  was  used. 

Figure  15  is  a  plot  of  the  total  moisture  in  the  top  five  feet  of 
soil  versus  time.  These  data  were  collected  using  the  neutron 
probe.  In  this  field  the  probe  readings  were  very  different  from 
each  other,  which  is  not  uncommon  for  this  method  of  soil  moisture 
monitoring . 

vi.   Upflux  and  salt  movement 

Upflux  is  defined  as  the  movement  of  water  from  the  saturated  zone 
into  the  crop  root  zone.  The  amount  of  upflux  which  occurs  in  a 
field  between  irrigations  can  be  estimated  by  subtracting  the 
observed  change  in  soil  moisture  from  the  estimated  ET .  Table  5 
summarizes  the  results  of  upflux  calculations  for  each  of  the  three 
neutron  probe  sites  in  field  P-1. 

The  average  upflux  based  on  Table  5  was  57%  of  ET.  This  means  that 
between  the  irrigations  57%  of  water  used  by  the  crop  was  pulled 
from  the  water  table.  Since  the  growing  season  ET  was  22.9  inches, 
13.1  inches  came  from  upflux. 
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Table  5.   Upflux  estimates  for  field  P-1 


Upflux       Upflux    Avg 
Site    Period        SM      ET    (ET-Change  SM)   %  of  ET   DTWT 


Change 

Period 

SM 

ET 

(in.) 

( in. 

5/26-6/9 

0.01 

0.68 

6/13-7/6 

1.54 

3.26 

7/11-7/22 

0.26 

3.02 

7/29-8/10 

0.48 

2.93 

8/12-9/29 

0.98 

6.91 

(%)    (ft.) 


5/26-6/9 

0.20 

0.68 

6/13-7/6 

2.34 

3.26 

7/11-7/22 

1.04 

3.02 

7/29-8/10 

1.81 

2.93 

8/15-9/29 

4.58 

6.48 

5/26-6/9 

0.05 

0.68 

6/13-7/6 

0.90 

2.29 

7/11-7/22 

2.14 

3.02 

7/25-8/10 

2.49 

4.17 

8/12-9/29 

2.48 

7.56 

OW-12   5/26-6/9     0.01    0.68        0.67  99 

1.72  53 

2.76  91 

2.45  84 

5.93  86 

Sum       3.27   16.80       13.53  81%      6.7 

OW-15   5/26-6/9     0.20    0.68        0.48  70 

0.92  28 

1.98  66 

1.12  38 

1.90  29 

Sum       9.97   16.37        6.40  39%      7.4 

OW-19   5/26-6/9     0.05    0.68        0.63  93 

1.39  61 

0.88  29 

1.69  40 

5.08  67 

Sum       8.47   17.72        9.66  52%     7.1 


It  is  important  to  note,  however,  that  over  the  course  of  the 
season,  deep  percolation  during  irrigations  was  greater  than  the 
upflux  (see  the  next  section).  Thus,  there  was  a  net  positive 
downward  movement  of  water.  The  total  upward  movement  of  salts 
from  the  water  table  based  upon  an  average  water  table  EC  of  6.58 
dS/m  calculates  to  635  tons.  IJ_  there  had  been  no  deep 
percolation,  the  average  ECe  of  the  top  five  feet  of  soil  would  be 
expected  to  increase  by  approximately  2  dS/m  due  to  upflux. 

According  to  the  theory  postulated  by  the  Phase  I  report,  the  net 
positive  downward  water  movement  also  insured  a  net  positive 
downward  salt  movement.  Thus,  it  is  expected  that  all  of  the  635 
tons  of  salt  which  moved  up  from  the  water  table  over  the  course 
of  a  year,  also  moved  back  down  into  the  water  table.  Root  zone 
soil  samples  will  be  taken  to  verify  or  dispute  this  contention. 
The  previous  calculation  showed  that  the  very  worst  short  term 
salinization  which  could  have  occurred  in  this  field  was  2  dS/m. 
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vii.   Water  destinations 

Figure  16  is  an  annual  water  destination  graph  for  Field  P-1.  All 
of  the  infiltrated  water  from  the  1988  crop  season  is  accounted  for 
in  this  graph.  Of  the  26.8  inches  of  infiltrated  irrigation  water 
plus  rain,  23.3  inches  evaporated  or  was  transpired  by  the  crop. 
The  net  deep  percolation  to  the  water  table  calculates  to  3.5 
inches.  Since  the  tile  drain  outflow  from  this  field  was  also  3.5 
inches,  the  net  deep  flow  was  0.0  inches.  This  means  that  over  the 
course  of  the  season,  no  more  water  entered  the  field  from  lateral 
or  deep  subsurface  inflow  than  left  the  field  through  lateral  or 
deep  subsurface  flow. 

The  amount  of  water  infiltrated  in  the  low  quarter  of  the  field 
was  estimated  at  25.1  inches.  Since  this  exceeded  ET  by  1.8 
inches,  there  was  a  net  positive  downward  movement  of  water  in  all 
of  the  field.  According  to  the  theory  postulated  in  the  Phase  I 
report,  the  goal  of  a  subsurface  drainage  reduction  program  should 
be  to  reduce  this  number  to  equal  the  theoretical  leaching 
requirement.  The  water  quality  of  Panoche  Water  District  Water  and 
the  crops  grown  in  the  field  give  a  theoretical  leaching 
requirement  of  less  than  5%  of  the  ET ,  or  1.2  inches.  This 
suggests  that  the  amount  of  drain  water  flowing  out  of  this  field 
in  1988  was  about  as  low  as  is  possible  without  creating  long  term 
salt  build  up  in  the  soil. 

viii.   DRAINMOD  simulations 

Figure  17  shows  graphs  of  the  depth  to  water  table  predicted  by 
DRAINMOD  simulations  of  tiles  number  4  and  5  compared  to  actual 
field  readings.  The  general  shape  of  the  DRAINMOD  curves  is 
similar  to  actual  readings  in  both  cases.  However,  DRAINMOD  under- 
predicted  the  build  up  in  the  water  table  immediately  following  the 
pre-irrigation  and  June  irrigation.  DRAINMOD  also  over-predicted 
the  water  table  response  to  May  rainfall  events.  The  actual 
readings  also  show  a  rise  in  the  water  table  caused  by  irrigations 
of  tomatoes  to  the  South  in  late  April  which  was  not  predicted  by 
DRAINMOD.  The  simulation  of  tile  number  5  showed  a  much  better 
correlation  than  did  the  simulation  of  tile  number  4. 

Figure  18  shows  graphs  of  the  tile  drain  outflow  predicted  by 
DRAINMOD  simulations  of  tiles  number  4  and  5  compared  to  actual 
field  readings.  The  general  shape  of  the  DRAINMOD  curves  is 
similar  to  actual  readings  in  both  cases.  However,  DRAINMOD  under- 
predicted  the  flow  rate  immediately  following  the  pre-irrigation 
and  June  irrigation.  DRAINMOD  also  over-predicted  the  tile  flow 
response  to  May  rainfall  events.  The  actual  readings  also  show  a 
flow  increase  caused  by  irrigations  of  tomatoes  to  the  South  in 
late  April  which  was  not  predicted  by  DRAINMOD.   Again,  the 
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simulation  of  tile  number  5  showed  a  much  better  correlation  than 
did  the  simulation  of  tile  number  4.  Table  6  compares  the  monthly 
outflows  of  the  DRAINMOD  simulations  versus  actual. 

The  inability  of  DRAINMOD  to  model  the  regional  groundwater  effects 
of  either  deep  downward  or  upward  flow  probably  accounts  for  much 
of  the  variation  between  predicted  and  actual  events.  Inaccuracies 
in  input  data  also,  no  doubt,  contribute  as  do  limitations  in  the 
model's  mathematics. 

Table  6.   Actual  Tile  Outflows  vs.  DRAINMOD,  Field  P-1 


Tile  5 

DRAINMOD 

Tile  4 

DRAINMOD 

Month 

Flow 

Flow 

Flow 

Flow 

(in) 

(in) 

(in) 

(in) 

Jan 

0.12 

0.33 

1.40 

0.00 

Feb 

4.47 

2.71 

1.56 

2.54 

Mar 

3.76 

3.85 

0.07 

3.13 

Apr 

1.86 

1.85 

0  .00 

0.94 

May 

0.56 

1.32 

0  .00 

1.13 

Jun 

2.06 

0.95 

0  .00 

0.92 

Jul  -  Dec 

0.00 

0.00 

0  .  00 

0.00 

Total 

12.83 

11.01 

3.03 

8.66 
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B.   Field  P-2 

i .   I rr igation 

Table  7  summarizes  data  collected  on  the  irrigations  of  Field  P-2 
in  1988. 

Table  7.   Field  P-2  Irrigation  Summary 


Start 

End 

Gross 

Lost 

Avg. 

Dist. 

Date 

Date 

Applied 

Runoff 

Inf ilt. 

Unif .* 

(dd/mm) 

(dd/mm) 

(in.) 

(in.) 

(in.) 

(%) 

02/09 

02/27 

12.29 

0.09 

12.2 

83 

06/09 

06/14 

5.32 

0.03 

5.3 

60 

07/04 

07/11 

7.78 

0.07 

7.7 

80 

07/22 

07/28 

6.14 

0.21 

5.9 

80# 

08/09 

08/14 

7.31 

0.25 

7.1 

80# 

Total 

for  year 

38.84 

0.65 

38.2 

78 

*  Average  distribution  uniformity  of  test  furrows 

#  Estimated  values,  as  they  were  not  measured 

This  field  was  deep  ripped  following  the  saf flower  harvest  in  1987. 
The  pre-season  irrigation  was  made  in  February.  As  a  result  of 
JML's  evaluation  of  the  performance  on  the  first  set,  the  farmer 
changed  from  48  hour  sets  with  low  furrow  flows  to  36  hour  sets 
with  higher  flows  for  the  remainder  of  the  irrigation.  This 
reduced  the  amount  of  infiltrated  water  by  an  unknown  amount.  The 
amount  which  infiltrated  was  still  much  more  than  was  necessary  to 
replenish  the  potential  root  zone  for  the  coming  season.  However, 
the  extremely  high  initial  soil  intake  characteristics  and  the  need 
to  adequately  "sub"  up  the  furrows  to  provide  seed  bed  moisture 
prevented  any  further  reductions. 

The  first  in-season  irrigation  was  in  June.  This  was  done  on  every 
other  furrow  to  reduce  the  amount  infiltrated.  Two  out  of  every 
three  rows  had  received  wheel  traffic,  while  the  remaining  was  a 
guess  row.  This  resulted  in  a  huge  difference  in  the  amount 
infiltrated  and  in  the  distribution  uniformity  (DU)  in  the  wheel 
rows  versus  the  guess  rows.  The  low  DU  and  relatively  high  amount 
infiltrated  for  the  field  was  due  to  the  guess  rows, 
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than  was  needed  to  replenish  soil  moisture  because  the  high  initial 
rate  prohibited  any  further  reductions. 

ii.   Tile  drain  flows,  quality,  and  loads 

Table  8  is  a  summary  of  the  tile  drain  data  for  the  1988  crop  year. 

The  flow  of  tile  water  and  the  constituent  loads  increase  from  the 
South  (number  1)  to  the  North  (number  2).  The  higher  tile  flows 
and  loads  came  from  the  areas  with  shallower  groundwater  levels. 
Figure  19  is  a  graph  of  tile  flows  versus  time.  The  tile  flows 
responded  strongly  to  all  of  the  irrigation  events.  They  also 
showed  a  response  in  March  to  the  irrigation  of  a  tomato  field  to 
the  south.  The  shapes  of  each  tile  drain  hydrograph  were  the  same 
except  that  the  amplitudes  varied,  increasing  from  South  to  North. 

Multiplying  the  water  analysis  results  for  Panoche  Water  District 
water  by  the  amount  of  infiltrated  water  gives  the  following  loads 
infiltrated  in  Field  P-2:  salt  =  283  tons,  boron  =  401  lbs, 
selenium  =  1.96  lbs. 

Comparing  the  incoming  to  outgoing  loads  shows  that  the  incoming 
salts  were  equal  to  68  percent  of  the  outgoing  salts.  Incoming 
boron  was  equal  to  20  percent  of  the  outgoing  boron.  Incoming 
selenium  was  equal  to  78  percent  of  the  outgoing  selenium. 

Table  8.   Field  P-2  Tile  Drain  Data  Summary 


Tile 

Total 

Avg. 

Avg. 

Avg. 

Salt 

Boron 

Se 

Number 

Flow 

EC 

Boron 

Se 

Load 

Load 

Load 

(#) 

(ac-ft) 

(dS/m) 

(mg/1) 

(mg/l) 

(tons) 

(lbs) 

(lbs) 

1 

6.2 

3.36 

4.60 

0.008 

19.0 

78.3 

0.13 

2 

10.3 

3.27 

5.15 

0.006 

31.8 

144.5 

0.16 

3 

17.4 

3.16 

5.34 

0.006 

51.7 

261.7 

0.27 

4 

26.8 

2.95 

4.26 

0.009 

73.8 

379.9 

0.57 

5 

27.6 

3.16 

5.48 

0.005 

81.7 

422.5 

0.41 

6 

25.8 

3.07 

5.61 

0.008 

73.4 

376.4 

0.46 

7* 

27.6 

3.16 

5.07 

0.007 

83.0 

380.6 

0.52 

Total 

142.9 

3.16 

5.07 

0.007 

414.4 

2043.9 

2.52 

(11.4  in.) 

*  For  this  analysis,  tile  7  is  the  deep  drain.  Actual  flows  and 
quality  could  not  be  measured  since  it  also  intercepts  seepage  from 
the  T  canal.  It  was  assumed  to  have  an  outflow  equal  to  tile  5  and 
EC,  boron,  and  selenium  equal  to  the  field  average. 
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iii.   Shallow  groundwater  levels  and  quality 

The  spatial  variation  of  the  water  table  is  illustrated  by  Figure 
20.  Figure  20  is  a  depth  to  groundwater  map  for  January  21,  1988. 
On  that  day,  the  depth  to  groundwater  varied  from  a  low  of  5.7  to 
a  high  of  greater  than  10  feet.  The  general  trend  of  decreasing 
depths  to  water  table  from  the  Southwest  to  the  Northeast  mean  that 
the  shallow  groundwater's  slope  is  not  as  steep  as  the  ground 
surface  slope  (see  Figure  21,  ground  surface  elevations).  The 
presence  of  the  deep  open  drain  along  the  eastern  field  edge  can 
be  seen  in  the  deeper  groundwater  contours  near  it. 
The  large  variation  in  the  water  table  across  this  field  was  one 
of  the  reasons  the  decision  was  made  to  postpone  controlled 
drainage  attempts  until  1989. 

The  time  variation  of  the  water  table  is  illustrated  in  Figure  22. 
This  is  a  plot  of  selected  observation  wells  during  the  1988  crop 
season.  Shown  is  the  significant  water  table  response  to  all  of 
the  irrigations  and  to  the  March  tomato  irrigation  on  the  field  to 
the  South. 

In  January  the  average  EC  was  3.16  ds/m  and  the  average  boron  was 
3.73  mg/1.  In  July  the  average  EC  was  4.15  dS/m  and  the  average 
boron  was  3.77  mg/1.  Figure  23  shows  the  spatial  variability  of 
the  groundwater  EC.  Figure  24  shows  the  spatial  variability  of  the 
groundwater  boron  content.  Both  show  a  very  similar  pattern. 
Values  tend  to  be  higher  in  the  east  half  where  the  water  table  is 
shallower  but  the  highest  values  do  not  coincide  with  the  areas  of 
shallowest  groundwater. 

iv.   Soils 

The  soils  in  the  top  ten  feet  of  the  field  were  found  to  be  pre- 
dominantly silty  clay-loam  and  silty  clay.  The  soils  tended  to  be 
less  variable  than  in  Field  P-1.  Some  layers  of  other  soil  types 
were  found,  from  loamy  sands  to  dense  sandy  clays.  Average  results 
of  the  soil  salinity  analysis  are  in  Table  9. 
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Table  9.   Average  Soil  EC  and  Boron,  Field  P-2 


Depth 

Avg  . 

Avg  . 

Range 

ECe 

Boron 

(ft) 

(dS/tn) 

(mg/l) 

0-1 

1.10 

0.95 

1-2 

1.05 

1.15 

2-3 

2.24 

2.70 

3-4 

3.04 

3.66 

4-5 

2.60 

3.52 

Avg  . 

2.00 

2.40 

The  soil  ECe  and  boron  varied  considerably  with  location  from  the 
above  averages.  None  of  the  sites  tested  had  average  ECe  or  boron 
levels  which  would  cause  yield  reductions  on  cotton.  Figure  25 
shows  lines  of  equal  average  root  zone  ECe  for  Field  P-2.  Figure 
26  shows  lines  of  equal  average  root  zone  boron  content  for  Field 
P-2.  The  ECe  and  soil  boron  patterns  were  nearly  identical.  The 
patterns  also  correlated  highly  with  the  depth  to  shallow 
groundwater  map.  Areas  with  shallower  groundwater  showed  higher 
ECe  and  soil  boron  while  areas  with  deeper  groundwater  showed  lower 
ECe  and  soil  boron. 

Figure  27  is  a  plot  based  upon  the  averages  of  all  the  soil 
moisture  characteristic  curves  derived  from  observations  wells  1, 
12,  13  &  20.  This  curve  is  very  similar  in  shape  to  Field  P-l's, 
but  is  shifted  higher,  perhaps  due  to  higher  clay  content  in  P-2's 
soils. 

In-field  tests  at  observation  wells  1,  12,  13,  and  20  gave 
saturated  horizontal  hydraulic  conductivity  values  of  1.1,  2.7, 
2.2,  and  4.1  in/hr  respectively. 

V.   Evapotranspiration  and  soil  moisture 

Figure  28  is  a  plot  of  the  evapotranspiration  potential  (ETP) ,  the 
ETP  multiplied  by  the  crop  coefficient  (ETP  x  Kc),  and  the  crop 
evapotranspiration  (CROP  ET) .  The  total  growing  season  CROP  ET  was 
21.1  inches.  The  total  growing  season  ETP  x  Kc  was  24.4  inches. 
Since  the  SWAP/ET  model  does  not  include  the  effect  of  upflux  from 
a  shallow  water  table,  the  actual  soil  dryness  was  less  than 
predicted  by  SWAP/ET.  For  this  reason,  actual  ET  numbers  should 
fall  somewhere  between  CROP  ET  and  ETP  x  Kc.  For  DRAINMOD 
simulations  on  this  field,  the  ETP  x  Kc  was  used  because  the  water 
table  was  relatively  shallow  during  the  whole  crop  season. 
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Figure  29  is  a  plot  of  the  total  moisture  versus  time  in  the  top 
five  feet  of  soil.  Calculations  of  total  crop  season  upflux  using 
the  estimated  ET  and  the  neutron  probe  readings  give  the  following 
numbers:  at  OWl ,  upflux  =  9.4  in.  (40%  of  ET) ;  at  0W12,  upflux  = 
14.4  in.  (62%  of  ET);  at  OW20,  upflux  =  7.7  in.  (33%  of  ET)  . 

vi.   Upflux  and  salt  movement 

Upflux  is  defined  as  the  movement  of  water  from  the  saturated  zone 
into  the  crop  root  zone.  The  amount  of  upflux  which  occurs  in  a 
field  between  irrigations  can  be  estimated  by  subtracting  the 
observed  change  in  soil  moisture  from  the  estimated  ET .  Table  10 
summarizes  the  results  of  upflux  calculations  for  each  of  the  three 
neutron  probe  sites  in  field  P-2. 

Table  10.   Upflux  estimates  for  field  P-2 


Change 

Upflux 

Upflux 

Avg  . 

Site 

Period 

SM 

ET 

(ET-Change  SM) 

%  of  ET 

DTWT 

(in.) 

(in.) 

(%) 

(ft.) 

OW-12 

5/26-6/17 

0.47 

1.14 

0.67 

59 

6/22-7/11 

1.50 

3.16 

1.66 

53 

7/13-7/22 

0.63 

2.50 

1.87 

75 

7/25-8/12 

2.11 

5.51 

3.40 

62 

8/15-9/29 

2.61 

8.73 

6.12 

70 

Sum 

7.32 

21.04 

13.72 

65% 

4.8 

OW-1 

5/26-6/17 

0.18 

1.14 

0.96 

84 

6/22-7/11 

2.22 

3.16 

0.94 

29 

7/13-7/22 

0.66 

1.90 

1.24 

65 

7/25-8/12 

2.80 

5.51 

2.71 

49 

8/15-9/29 

6.27 

8.73 

2.46 

28 

Sum 

12.14 

20.44 

8.31 

41% 

6.6 

OW-20 

5/26-6/17 

0.91 

1.14 

0.23 

20 

6/22-7/11 

4.23 

3.16 

-1.07 

-34 

7/13-7/22 

1.37 

2.50 

1.13 

45 

7/25-8/12 

3.07 

5.51 

2.44 

44 

8/15-9/29 

4.58 

8.73 

4.15 

48 

Sum 

14.16 

21.04 

6.88 

33% 

3.8 

The  average  upflux  based  on  Table  10  was  46%  of  ET .  This  means 
that  between  the  irrigations  46%  of  water  used  by  the  crop  was 
pulled  from  the  water  table.  Since  the  growing  season  ET  was  22.9 
inches,  10.5  inches  came  from  upflux. 
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It  is  important  to  note,  however,  that  over  the  course  of  the 
season,  deep  percolation  during  irrigations  was  greater  than  the 
upf lux  (see  the  next  section) .  Thus,  there  was  a  net  positive 
downward  movement  of  water.  The  total  upward  movement  of  salts 
from  the  water  table  based  upon  an  average  water  table  EC  of  4.15 
dS/m  calculates  to  563  tons.  Ij_  there  had  been  no  deep 
percolation,  the  average  ECe  of  the  top  five  feet  of  soil  would  be 
expected  to  increase  by  approximately  1.1  dS/m  due  to  upf lux. 

According  to  the  theory  postulated  by  the  Phase  I  report,  the  net 
positive  downward  water  movement  also  insured  a  net  positive 
downward  salt  movement.  Thus,  it  is  expected  that  all  of  the  563 
tons  of  salt  which  moved  up  from  the  water  table  over  the  course 
of  a  year,  also  moved  back  down  into  the  water  table.  Root  zone 
soil  samples  will  be  taken  to  verify  or  dispute  this  contention. 
The  previous  calculation  showed  that  the  very  worst  short  term 
salinization  which  could  have  occurred  in  this  field  was  1.1  dS/m. 

vii.   Water  destinations 

Figure  30  is  an  annual  water  destination  graph  for  Field  P-2.  This 
graph  accounts  for  all  of  the  infiltrated  water  for  the  1988  crop 
season.  Of  the  43  inches  of  infiltrated  irrigation  water  plus 
rain,  24.8  inches  evaporated  or  was  transpired  by  the  crop.  The 
deep  percolation  calculates  to  18.2  inches.  Since  the  tile  drain 
outflow  from  this  field  was  11.2  inches,  the  net  subsurface  outflow 
other  than  tile  flow  was  7.0  inches. 

The  amount  of  water  infiltrated  in  the  low  quarter  of  the  field 
was  estimated  at  34.6  inches.  Since  this  exceeded  ET  by  9.8 
inches,  there  was  a  net  positive  downward  movement  of  water  in  all 
of  the  field.  The  9.8  inch  figure  is  much  more  than  the 
theoretical  leaching  requirement  of  approximately  1.2  inches.  This 
means  that  8.6  inches  of  subsurface  drainage  could  have  been 
reduced  without  causing  a  long  term  salt  buildup  in  this  field. 

viii.   DRAINMOD  simulations 

Figure  31  shows  graphs  of  the  depth  to  water  table  (for  field  P-2) 
predicted  by  DRAINMOD  simulations  of  tiles  number  2  and  4,  compared 
to  the  field  readings.  The  general  shape  of  the  DRAINMOD  curves 
is  similar  to  the  field  readings  in  both  cases.  However,  DRAINMOD 
under-predicted  the  build  up  in  the  water  table  immediately 
following  the  irrigations.  DRAINMOD  also  over-predicted  the  water 
table's  response  to  May  rainfall  events.  The  field  data  also  show 
a  rise  in  the  water  table  caused  by  irrigations  of  tomatoes  to  the 
South  in  late  April.  This  was  not  predicted  by  DRAINMOD.  The 
simulation  of  tile  number  4  showed  a  much  better  correlation  than 
did  the  simulation  of  tile  number  2. 
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Figure  32  shows  graphs  of  the  tile  drain  outflow  predicted  by 
DRAINMOD  simulations,  for  tiles  number  2  and  4,  compared  to  the 
field  readings.  The  general  shape  of  the  DRAINMOD  curves  is 
similar  to  the  field  data  in  both  cases.  However,  DRAINMOD  under- 
predicted  the  flow  rate  immediately  following  the  irrigations. 
DRAINMOD  also  over-predicted  the  tile  flow  response  to  May  rainfall 
events.  The  actual  readings  also  show  a  flow  increase  caused  by 
irrigations  of  tomatoes  to  the  South  in  late  April,  which  was  not 
predicted  by  DRAINMOD.  Again,  the  simulation  of  tile  number  4 
showed  a  much  better  correlation  than  did  the  simulation  of  tile 
number  2.  Table  11  compares  the  monthly  outflows  of  the  DRAINMOD 
simulations  versus  actual. 

Table  11.   Actual  Tile  Outflows  vs.  DRAINMOD,  Field  P-2 


Tile  4 

DRAINMOD 

Tile  2 

DRAINMOD 

Month 

Flow 

Flow 

Flow 

Flow 

(in) 

(in) 

(in) 

(in) 

Jan 

0.18 

0.11 

0.03 

0.00 

Feb 

2.58 

1.72 

1.06 

1.31 

Mar 

1.65 

2.05 

0.71 

1.66 

Apr 

1.15 

0.87 

0.38 

0.71 

May 

0.76 

0.53 

0.15 

0.47 

Jun 

1.42 

1.00 

0.37 

0.97 

Jul 

3.56 

3.86 

2.07 

3.79 

Aug 

2.16 

2.65 

0.54 

2.66 

Sep 

0.36 

0.02 

0.00 

0.03 

Total 

13.82 

12.81 

5.31 

11.60 

The  inability  of  DRAINMOD  to  model  the  regional  groundwater  effects 
of  either  deep  downward  or  upward  flow  probably  accounts  for  much 
of  the  variation  between  predicted  and  actual  events.  Inaccuracies 
in  input  data  also  no  doubt  contribute,  as  do  limitations  in  the 
model's  mathematics. 
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C.   Field  B10-2 

i.   Irrigation 

Table  12  summarizes  data  collected  on  the  irrigations  of  Field  B10- 
2  in  1988. 

Table  12.   Field  B10-2  Irrigation  Summary 


Start 
Date 
(dd/mm) 

End 
Date 
(dd/mm) 

Gross 

Applied 

(in.) 

Run- 
off 
(in.) 

Avg  . 
Inf ilt. 
(in.) 

Dist. 
Unif .* 
(%) 

04/04 
06/16 
07/02 
07/22 
08/09 

04/17 
06/20 
07/08 
07/28 
08/14 

12.5 
4.1 
1.4 
1.2 
0.5 

4.1 
1.4 
0.6 
0.7 
0.3 

8.4 
2.7 
0.8 
0.6 
0.2 

87 

87 

92 

92# 

92# 

Total 

for  year 

19.7 

7.1 

12.7 

88 

*  Average  distribution  uniformity  of  test  furrows 

#  Estimated  values,  as  they  were  not  measured 

The  pre-season  irrigation  was  made  in  April.  An  evaluation  of  the 
performance  on  the  first  set  found  that  changing  from  48  hour  sets 
with  low  furrow  flows  to  24  hour  sets  with  higher  flows  for  the 
remainder  of  the  irrigation  would  have  reduced  the  deep  percolation 
by  1.3  inches.  Since  additional  lost  runoff  made  the  practice  no 
more  economical  for  the  farm,  the  farmer  stayed  with  the  48  hour 
sets.  The  amount  which  infiltrated  was  much  more  than  was 
necessary  to  replenish  the  potential  root  zone  for  the  coming 
season.  However,  the  soil  intake  characteristics  and  the  need  to 
adequately  "sub"  up  the  furrows  to  provide  seed  bed  moisture 
prevented  major  reductions  in  deep  percolation. 

All  of  the  in-season  irrigations  were  made  with  two  passes  across 
the  field  using  wheel  compacted  furrows  only  (one  out  of  every 
three  furrows  was  a  non-wheel  row  and  was  not  irrigated).  This 
not  only  reduced  the  amount  of  deep  percolation,  but  also  helped 
to  prevent  crop  damage  caused  by  saturation  of  the  melon  beds.  In 
previous  years,  melons  in  this  field  had  suffered  serious  crop 
losses  due  to  this  problem.  The  crop  yields  this  year  were  much 
higher  than  in  previous  years.  During  the  last  three  irrigations 
the  soil  had  sealed  up  enough  to  where  the  soil  moisture  was  not 
being  completely  replenished. 
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ii.   Tile  drain  flows,  quality,  and  loads 

Figure  33  shows  the  time  variation  in  tile  system  flows  for  Field 
B10-2.  Table  13  is  a  summary  of  the  tile  drain  data  for  the  1988 
crop  year. 

Table  13.   Field  B10-2  Tile  Drain  Data  Summary 


Field 
Number 

(#) 

Total 
Flow 
(ac-ft) 

Avg. 
EC 
(dS/m) 

Avg  . 
Boron 
(mg/l) 

Avg  . 
Se 
(mg/l) 

0.102 
0.054 

0.074 

Salt 

Load 

(tons) 

324 
418 

741 

Boron 
Load 

(lbs) 

Se 
Load 

(lbs) 

B10-1 
B10-2 

42.0 
57.9* 

99.9 

7.88 
7.37 

7.57 

5.03 
6.29 

5.76 

575 
990 

11.7 
8.5 

Total 

1565 

20.2 

4.8  inches 


Average  laboratory  results  for  Broadview  Water  District's 
irrigation  water,  weighted  by  the  amount  infiltrated  during  each 
irrigation  were:  EC  =  0.85  dS/m  (544  mg/l),  boron  =  0.62  mg/l,  and 
se  =  0.0054  mg/l.  The  calculated  load  of  constituents  in  the 
infiltrated  water  was,  therefore:  salt  =  114  tons,  boron  =  258  lbs, 
selenium  =  2.2  lbs. 

Comparing  the  incoming  to  outgoing  loads  for  Field  B10-2  shows  that 
the  outgoing  salts  were  3.7  times  higher  than  the  incoming  salts, 
outgoing  boron  was  3.8  times  higher  than  incoming  boron,  and 
outgoing  selenium  was  3.9  times  higher  than  incoming  selenium. 

iii.   Shallow  groundwater  levels  and  quality 

The  spatial  variation  of  the  water  table  is  illustrated  by  Figure 
34.  This  figure  is  a  depth  to  groundwater  map  for  March  15,  1988, 
which  was  prior  to  the  pre-irrigation.  On  that  day,  the  depth  to 
groundwater  varied  from  a  low  of  5.9  to  a  high  of  8.2  feet.  The 
groundwater  tended  to  get  deeper  from  the  South  to  the  North.  This 
was  possibly  due  to  the  Broadview  Water  District  irrigation  ditch 
which  runs  along  the  south  edge  of  the  field.  Figure  35  shows  the 
ground  surface  elevations  of  the  field. 

The  time  variation  of  the  water  table  is  illustrated  in  Figure  36. 
This  is  a  plot  of  selected  observation  wells  through  the  1988  crop 
season.  Shown  is  how  the  water  table  responded  significantly  to 
the  pre-season  irrigation  and  first  in-season  irrigation  but  had 
little  response  to  the  remaining  irrigations.  The  water  table 
showed  little  observable  response  to  the  irrigation  of  any  of  the 
surrounding  fields  (which  are  all  tiled). 
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In  February,  the  average  EC  was  8.49  dS/m  and  the  average  boron 
was  6.33  mg/1.  Figures  37  and  38  show  the  spatial  variability  of 
the  groundwater  EC  and  boron  content,  respectively.  The  patterns 
were  very  similar  except  for  the  southeast  quarter  of  the  field. 
No  correlation  with  water  table  depths  or  soil  types  could  be  found 
to  explain  the  patterns. 

iv.   Soils 

The  soils  in  the  top  ten  feet  of  the  field  were  found  to  be 
stratified  and  vary  considerably  from  hole  to  hole.  The  dominant 
soil  type  was  silty  clay-loam.  Soil  types  observed  in  layers  in 
many  of  the  holes  ranged  from  loamy  sand  to  dense  clay.  In  general 
the  soils  tended  to  become  finer  from  the  Northwest  to  the 
Southeast.  The  northwest  quarter  of  the  field  has  sandy  layers  in 
most  of  the  profiles.  Average  results  of  the  soil  salinity 
analysis  are  in  Table  14. 

Table  14.   Average  Soil  EC  and  Boron,  Field  B10-2 

Avg . 
Boron 
(mg/l) 


Def 

)th 

Rar 

iqe 

(ft) 

0- 

■1 

1- 

•2 

2- 

■3 

3- 

•4 

4- 

•5 

Av 

"3- 

EC 

.e 

(dS/m) 

2. 

55 

2, 

50 

3. 

39 

3, 

89 

4. 

34 

1 

.15 

2 

.15 

3 

.13 

3 

.03 

2 

.82 

Avg.        3.3  3        2.4  6 


The  soil  ECe  and  boron  varied  considerably  with  location  from  the 
above  averages.  Some  of  the  sites  tested  had  average  ECe  or  boron 
levels  which  could  cause  yield  reductions  on  melons.  Figure  39 
shows  lines  of  equal  average  root  zone  ECe  for  Field  B10-2.  Figure 
40  shows  lines  of  equal  average  root  zone  boron  content. 

In-field  tests  at  observation  wells  1,  7,  16  and  20  gave  saturated 
horizontal  hydraulic  conductivity  values  of  1.2,  1.8,  0.8  and  0.5 
in/hr  respectively. 

V.   Evapotranspiration  and  soil  moisture 

Figure  41  is  a  plot  of  the  evapotranspiration  potential  (ETP)  ,  the 
ETP  multiplied  by  the  crop  coefficient  (ETP  x  Kc)  ,  and  the  crop 
evapotranspiration  (CROP  ET) .  The  effect  of  deficit  irrigations 
with  no  consideration  of  upflux  make  the  CROP  ET  figures  from 
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Figure  42  is  a  plot  of  the  total  moisture  in  the  top  five  feet  of 
soil  versus  time.  In  this  field  the  probe  readings  probably  did 
not  measure  all  of  the  soil  moisture  changes  which  were  taking 
place  in  the  field.  This  is  because  the  access  tubes  were  placed 
in  the  center  of  the  80  inch  wide  melon  bed  and  could  not  record 
moisture  changes  near  the  furrows.  Three  tubes  per  bed,  with  two 
tubes  near  the  shoulder  and  one  in  the  middle,  would  be  a  better 
arrangement . 

vi .   Upflux  and  salt  movement 

Upflux  is  defined  as  the  movement  of  water  from  the  saturated  zone 
into  the  crop  root  zone.  The  amount  of  upflux  which  occurs  in  a 
field  between  irrigations  can  be  estimated  by  subtracting  the 
observed  change  in  soil  moisture  from  the  estimated  ET .  Table  15 
summarizes  the  results  of  upflux  calculations  for  each  of  the  three 
neutron  probe  sites  in  field  B10-2. 

The  average  upflux  based  on  Table  15  was  55%  of  ET.  This  means 
that  between  the  irrigations  55%  of  water  used  by  the  crop  was 
pulled  from  the  water  table.  Since  the  growing  season  ET  was  22.9 
inches,  12.6  inches  came  from  upflux. 
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Table  15.   Upflux  estimates  for  field  B10-2 

Change  Upflux  Upflux  Avg . 

Site    Period        SM      ET    (ET-Change  SM)  %  of  ET  DTWT 

(in.)   (in.)  (%)  (ft.) 


OW-1    5/25-6/9  0.70  0.645  0.055  0 

6/17-7/4  0.00  1.31  1.31  100 

7/8-7/25  1.42  4.55  3.13  69 

7/27-8/12  0.73  1.42  0.69  49 


Sum       2.85    7.93        5.08  64%      6.3 

OW-7    5/25-6/9     1.07    1.22        0.15  12 

1.56  65 

3.49  77 

0.42  30 


5/25-6/9 

1.07 

1.22 

6/17-7/4 

0.83 

2.39 

7/8-7/25 

1.06 

4.55 

7/27-8/12 

1.00 

1.42 

Sum       3.96    9.58        5.62  59%      5.7 

OW-20   5/25-6/9      .62    1.22        0.60  49 

0.31  13 

3.01  66 

0.14  10 


5/25-6/9 

.62 

1.22 

6/17-7/4 

2.08 

2.39 

7/8-7/25 

1.54 

4.55 

7/27-8/12 

1.28 

1.42 

Sum       5.52    9.58        4.06  42%     5.2 


It  is  important  to  note,  however,  that  over  the  course  of  the 
season,  deep  percolation  during  irrigations  was  greater  than  the 
upflux  (see  the  next  section).  Thus,  there  was  a  net  positive 
downward  movement  of  water.  The  total  upward  movement  of  salts 
from  the  water  table  based  upon  an  average  water  table  EC  of  8.49 
dS/m  calculates  to  743  tons.  IJ_  there  had  been  no  deep 
percolation,  the  average  ECe  of  the  top  five  feet  of  soil  would  be 
expected  to  increase  by  approximately  1.5  dS/m  due  to  upflux. 

According  to  the  theory  postulated  by  the  Phase  I  report,  the  net 
positive  downward  water  movement  also  insured  a  net  positive 
downward  salt  movement.  Thus,  it  is  expected  that  all  of  the  743 
tons  of  salt  which  moved  up  from  the  water  table  over  the  course 
of  a  year,  also  moved  back  down  into  the  water  table.  Root  zone 
soil  samples  will  be  taken  to  verify  or  dispute  this  contention. 
The  previous  calculation  showed  that  the  very  worst  short  term 
salinization  which  could  have  occurred  in  this  field  was  1.5  dS/m. 
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vii.  Water  destinations 

Figure  43  is  an  annual  water  destination  graph  for  Field  B10-2 
based  upon  the  data  collected.  This  graph  accounts  for  all  of  the 
infiltrated  water  from  the  1988  crop  season.  Of  the  16.1  inches 
of  infiltrated  irrigation  water  plus  rain,  12.2  inches  evaporated 
or  was  transpired  by  the  crop.  The  net  deep  percolation  to  the 
water  table  calculates  to  3.9  inches.  Since  the  tile  drain  outflow 
from  this  field  was  4.8  inches,  the  net  deep  flow  was  -0.9  inches. 
This  suggests  that  over  the  course  of  the  season,  slightly  more 
water  entered  the  field  from  lateral  or  deep  subsurface  inflow  than 
left  the  field  through  lateral  or  deep  subsurface  flow. 

The  amount  of  water  infiltrated  in  the  low  quarter  of  the  field 
was  estimated  at  14.5  inches.  Since  this  exceeded  ET  by  2.3 
inches,  there  was  a  net  positive  downward  movement  of  water  in  all 
of  the  field.  Assuming  a  leaching  requirement  of  less  than  5%  of 
the  ET  ,  or  0.6  inches,  the  amount  of  drainage  flowing  out  of  this 
field  in  1988  could  have  been  reduced  an  additional  1.7  inches 
without  creating  long  term  salt  buildup  in  the  soils. 

viii.   DRAINMOD  simulations 

Figure  44  shows  a  graph  of  the  depth  to  water  table  predicted  by 
DRAINMOD  simulations  compared  to  the  field  readings.  The  shape  of 
the  DRAINMOD  curve  is  very  similar  to  the  field  data.  However,  the 
DRAINMOD  predicted  water  table  drifted  increasingly  below  the 
actual  water  table  through  the  crop  season.  This  could  be  due  to 
overestimates  of  ET ,  or,  to  upward  flows  from  below  the  tile  lines 
which  are  not  adequately  modeled  in  DRAINMOD. 

Figure  45  shows  a  graph  of  the  tile  drain  outflow  predicted  by  the 
DRAINMOD  simulation  versus  recorded  tile  flow.  The  DRAINMOD  curve 
does  not  match  very  well,  especially  late  in  the  season.  The  flows 
fluctuated  more  following  pre-irrigation  than  those  represented  by 
DRAINMOD.  Some  of  this  can  be  explained  by  the  method  of  obtaining 
total  drain  flow.  This  is  done  by  subtracting  estimated  flows  from 
Field  B10-1  (based  on  instantaneous  measurements)  from  the 
totalized  sump  flows  to  obtain  Field  B10-2  flows. 
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Table  16.   Actual  Tile  Outflows  vs.  DRAINMOD,  Field  B10-2 


Actual 

DRAINMOD 

Month 

Flow 

Flow 

(in) 

(in) 

Jan 

0.00 

0.00 

Feb 

0.00 

0.00 

Mar 

0.00 

0  .00 

Apr 

2.18 

2.34 

May 

1.05 

1.66 

Jun 

1.16 

0.46 

Jul 

0.47 

0.00 

Aug-Dec 

0.00 

0.00 

Total 

4.86 

4.46 

Again,  the  inability  of  DRAINMOD  to  model  the  regional  groundwater 
effects  of  either  deep  downward  or  upward  flow  probably  accounts 
for  much  of  the  variation  between  predicted  and  actual  events. 
Inaccuracies  in  input  data  also  no  doubt  contribute,  as  do 
limitations  in  the  model's  mathematics. 
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IV.   CONCLUSIONS 


A.  Field  P-1 

The  combination  of  irrigation  management  and  controlled  drainage 
in  Field  P-1  resulted  in  the  total  tile  drain  outflows  being  very 
near  the  target  level  for  reduced  drainage  while  maintaining  a  long 
term  salt  balance. 

There  was  a  great  difference  between  the  outflows  of  the  individual 
tile  drains.  This  could  be  attributed  mainly  to  differences  in  the 
placement  depth.  For  the  shallow  tiles,  the  deep  percolation  from 
rains  and  irrigation  was  greater  than  the  tile  outflow.  For  the 
deeper  tiles,  the  deep  percolation  was  not  as  great  as  the  tile 
outflow.  On  the  field  as  a  whole,  the  deep  percolation  was 
approximately  equal  to  the  tile  outflow. 

Drainmod  simulations  were  successful  in  predicting  the  major  trends 
of  water  table  and  tile  outflow.  However,  the  effects  of  the 
irrigation  of  neighboring  fields  and  both  downward  and  upward  flow 
below  the  tile  drains  made  for  significant  differences  between 
predicted  and  actual  results.  With  further  refinement  of  the  input 
data  and  some  program  modifications,  DRAINMOD  should  be  a  useful 
tool  for  planning  the  drainage  reduction  program  in  this  field. 

Further  refinements  on  the  drainage  control  program  are  needed  on 
this  field.  Another  line  of  tile  valves  is  needed  to  assure  more 
uniform  water  table  control.  The  weir  settings  should  be  set  to 
have  the  same  depth  below  ground  surface  (rather  than  depth  above 
drain  invert  as  was  tried  this  year)  throughout  the  field.  This 
would  be  in  order  to  offset  the  effect  of  the  variations  in  drain 
depths.  The  weirs  also  need  to  be  set  higher  during  the  early  part 
of  the  year,  when  almost  all  of  the  deep  percolation  occurs. 

B.  Field  P-2 

The  irrigation  management  program  implemented  this  year  in  Field 
P-2  was  not  sufficient  to  reduce  tile  drain  flows  down  to  the 
minimum  level  allowed  for  long  term  salt  balance.  The  high  water 
intake  rates  throughout  the  year,  especially  in  furrows  without 
wheel  traffic,  were  primarily  responsible  for  excessive  deep 
percolation. 

Tile  drain  flows  varied  considerably  from  tile  to  tile.  This  could 
be  explained  by  a  variation  in  the  depth  to  the  water  table  across 
the  field.  The  water  table  variation  appears  to  be  caused  by  the 
ground  surface  slope  being  greater  than  the  water  table  slope. 
Tile  drains  in  the  deeper  water  table  area  had  outflows  much  less 
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than  the  deep  percolation.  Tile  drains  in  the  shallower  water 
table  area  had  outflows  closer  to  the  deep  percolation. 

DRAINMOD  simulations  were  successful  in  predicting  the  general 
trend  of  the  water  table  and  drain  flows.  However,  the  effect  of 
irrigations  in  neighboring  fields,  and  the  flow  of  subsurface  water 
below  the  tile  drains,  made  for  significant  differences  between 
predicted  and  actual  water  table  levels  and  flows. 

The  key  to  a  successful  water  table  management  program  in  this 
field  will  be  to  handle  the  relatively  large  variation  in  the  water 
table  depth.  This  may  require  splitting  the  field  into  perhaps  as 
many  as  four  different  units,  each  with  different  management 
schemes. 


C.   Field  B10-2 

The  irrigation  management  of  Field  B10-2  resulted  in  reductions  in 
deep  percolation.  However,  the  tile  drain  flows  were  still  higher 
than  the  target  level  for  drainage  reduction  with  long  term  salt 
balance. 

The  DRAINMOD  simulations  correlated  well.  However,  DRAINMOD 
increasingly  tended  to  overestimate  the  depth  to  the  water  table 
as  the  season  progressed.  Either  estimates  of  ET  were  too  high  or 
there  is  an  upward  flow  from  beneath  the  tile  drains. 

Successful  water  table  control  in  this  field  will  be  enhanced  by 
being  able  to  calibrate  DRAINMOD  better  to  model  more  closely  the 
behavior  of  the  water  table.  Since  this  field  had  less  spacial 
variability  in  the  water  table,  and  less  influence  from  surrounding 
fields,  it  should  be  the  easiest  of  the  three  fields  for  which  to 
accomplish  water  table  management. 


D.   General 

All  of  the  data  collected  under  this  phase  of  our  study  confirms 
and  supports  the  concept  put  forward  under  Phase  I  that  irrigation 
management  and  controlled  drainage  techniques  can  be  used  to  reduce 
the  flow  and  constituent  loads  from  subsurface  drainage  systems. 
JMLord's  position  concerning  the  feasibility,  applicability  and 
limitations  of  irrigation  management/controlled  drainage  techniques 
for  subsurface  drainage  reduction  has  been  written  up  in  detail  in 
the  Phase  I  report  to  the  S.J.V.D.P.  Results  of  field  trials  have 
reinforced  that  perspective.  It  has  become  clear  that  the 
application  of  the  techniques  will  require  developing  management 
plans  on  an  individual  field  basis;  and  that  an  understanding  of 
each  field's  dynamics  must  be  developed  before  the  plan  can  be 
made.   The  degree  of  variation  in  soils,  water  table,  tile  drain 
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configurations,  etc.  within  the  field  will  determine  the  practical 
limit  to  which  drainage  can  be  reduced. 

The  degree  of  correlation  which  was  achieved  between  DRAINMOD 
simulations  and  the  field  data  is  encouraging.  With  proper  input 
data,  this  study  has  shown  that  the  current  version  of  DRAINMOD 
can  do  an  adequate  job  of  describing  the  major  hydrologic  processes 
which  occur  in  a  tile  drained  field,  if  it  is  isolated  from 
surrounding  influences.  The  model,  therefore,  holds  promise  for 
use  as  a  tool  in  evaluating  the  effects  of  management  changes 
within  the  field.  However,  the  inability  of  the  current  version 
to  model  regional  effects  appears  to  be  the  biggest  limitation  to 
its  accuracy.  The  magnitude  of  the  regional  effects  is  very  site 
specific.  In  some  cases  the  regional  effects  are  small  enough  so 
that  the  current  version  of  DRAINMOD  is  adequate.  In  other  cases 
the  regional  effects  may  be  strong  enough  to  make  it  inadequate. 

A  new  version  of  DRAINMOD  (to  be  released  in  1989)  will  have  the 
capability  of  modeling  deep  downward  or  upward  flow  as  well  as 
lateral  subsurface  flows.  It  remains  to  be  seen  whether  the  new 
version  will  be  useful  for  modeling  the  regional  effects. 

DRAINMOD  has  been  found  to  show  promise  as  a  tool  for  developing 
the  management  plans.  However,  there  are  numerous  programming 
refinements  needed  to  make  it  easier  to  use  and  improve  its  ability 
to  simulate  field  conditions.  The  following  modifications  are 
suggested: 

1.  A  more  user-friendly  interface  is  needed  which  would 
allow  direct  entry  of  an  irrigation  schedule  and  rainfall 
events,  direct  importation  of  potential  ET  and  crop  curve 
files,  direct  transfer  of  data  from  the  soil  preparation 
programs  (which  take  soil  data  and  develop  the  necessary 
inputs  for  DRAINMOD) ,  and  the  creation  of  output  files 
in  the  Lotus  format  for  easy  graphing  and  comparison. 

2.  DRAINMOD' s  program  needs  to  be  modified  to  better  model 
the  moisture  status  of  the  crop  root  zone  and  upflux, 
calculate  the  soil  dryness  coefficient  for  reducing  crop 
ET  within  DRAINMOD  itself,  model  deep  upward  or  downward 
flows  below  the  tile  drains,  and  model  salt  flows  and 
loads-. 

E.   Applications  and  costs  of  irrigation  management/controlled 
drainage 

Irrigation  management/controlled  drainage  techniques  are  an 
emerging  technology.  Work  is  in  progress  toward  making  the 
techniques  into  a  service  which  can  be  offered  to  growers.  Exactly 
how  the  techniques  will  be  applied  and  at  what  cost  have  not  been 
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fully  determined  yet. 

Hardware  costs  are  the  easiest  to  determine.  For  the  three  fields 
studied,  the  costs  for  installing  weir  control  valves  are  $73,  $67 
and  $47  per  acre.  Costs  will  be  site  specific.  Steeper  slopes  and 
narrower  tile  spacings  will  have  the  highest  costs.  Flow  control 
valves  have  not  been  developed  yet  for  tile  drains.  The  hardware 
cost  of  flow  control  is  therefore  not  known. 

Consulting  costs  are  more  difficult  to  determine.  JMLord  has  a 
goal  of  getting  costs  of  annual  consulting  fees  down  to  $2,500.00 
per  field.  This  would  not  include  costs  for  installing  flow 
measurement  devices  for  irrigation,  tile  drain  or  tailwater  runoff. 
Nor  would  it  include  laboratory  tests  for  soil  or  water  analysis. 
These  costs  are  not  included  because  they  are  likely  to  be  required 
with  or  without  the  consulting  service.  In  many  cases  they  will 
be  provided  by  the  irrigation/water  or  drainage  district.  Where 
not  covered  by  a  district,  these  could  be  provided  by  the 
consultant  at  additional  cost. 
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V.  1989  FIELD  STUDIES 


A.   Field  P-1 

Field  P-1  will  be  in  cotton  again  in  1989.  Another  line  of  tile 
drain  valves  will  be  installed.  The  weirs  will  be  set  based  upon 
farther  DRAINMOD  "what  if"  scenarios.  Monitoring  will  be  done  as 
it  was  in  1988. 


B.   Field  P-2 

Field  P-2  will  be  farmed  with  75  acres  of  cotton  and  75  acres  of 
tomatoes  in  1989.  The  split  will  be  made  on  an  East-West  line. 
Because  the  cropping  split  and  the  extreme  water  table  variations 
would  make  water  table  control  very  difficult,  a  different  field 
was  selected.  The  new  field  is  a  150  acre  cotton  field  in 
Westlands  Water  District.  Controlled  drainage  will  be  implemented 
in  1989. 


C.  Field  B10-2 

Field  B10-2  will  be  farmed  in  tomatoes  in  1989.  Since  Field 
B10-1  was  in  tomatoes  in  1988,  irrigation  and  drainage  data  from 
it  will  be  used  to  develop  water  table  control  plans.  The  half 
mile  field  will  be  split  into  three  880  foot  furrow  runs.  Tile 
drain  valves  will  be  put  into  the  tile  laterals  two-thirds  of  the 
way  from  the  head  end  and  the  weirs  set  for  water  table  control. 

D.  DWR  Study 

The  California  Department  of  Water  Resources  has  recently  awarded 
a  contract  to  JMLord,  Inc.  for  a  study  on  shallow  groundwater 
management.  The  primary  difference  between  the  DWR  study  and  this 
one  is  that  the  DWR  study  will  focus  mostly  upon  developing  and 
testing  new  tile  drainage  systems  to  incorporate  shallow 
groundwater  management  capabilities.  This  study  focuses  on 
retrofitting  existing  systems  for  shallow  groundwater  management. 
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APPENDICES 

Appendix  A  Tile  Drain  Valve  Literature 

specifications  and  costs  of  Dos-Ir  Valves 

Appendix  B  Tile  Drain  Flows  and  Loads 

-  tile  water  flows 

-  laboratory  results  for  boron, 
selenium,  EC,  salt 

Appendix  C  Depth  to  Groundwater  Reading 

-  dates  and  observation  well  readings 

Appendix  D  Shallow  Groundwater  Lab  Analysis 

-  Observation  well  water  samples,  analyzed 
for  EC  and  boron 

Appendix  E  Soil  Profile  Logs 

-  Logs  of  soil  texture,  apparent  permeability, 
soil  color  and  remarks 

Appendix  F  Soil  Salinity  Analysis 

-  Analysis  of  five  foot  soil  profile  for  pH , 
ECe,  B,  Ca+Mg,  Na  and  SAR 

Appendix  G  Soil  Mositure  Characteristic  Curve  Data 

-  Results  of  pressure  plate  analysis 

Appendix  H  Neutron  Probe  Readings 

-  soil  mositure  monitoring  results 

Appendix  I  DRAINMOD  Input  Data 

-  example  of  input  data  required  to  run 
DRAINMOD 


Appendices  are  available  from: 
San  Joaquin  Valley  Drainage  Program 
2880  Cottage  Way,  Room  W-2143 
Sacramento,  CA  95825-1898. 
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Figure   17a.     Field  P-1    Comparison 
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Figure  17b.     Reld  P-1    Comparison 
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Figure  19b.     Field  P-2 
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